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These notes are prepared for a set of lectures delivered at the The Jfii CTP 
Workshop on Statistical Physics: "Dynamics of Fluctuating Interfaces and Re- 
lated Phenomena" , at Seoul National University, Korea. The lectures examine 
several problems related to non-equilibrium fluctuations of interfaces and flux 
lines. The first two introduce the phenomenology of depinning, with partic- 
ular emphasis on interfaces and contact lines. The role of the anisotropy of 
the medium in producing different universality classes is elucidated. The last 
two lectures focus on the dynamics of lines, where transverse fluctuations are 
also important. We shall demonstrate how various non-linearities appear in 
the dynamics of driven flux lines. The universality classes of depinning, and 
also dynamic roughening, are illustrated in the contexts of moving flux lines, 
advancing crack fronts, and drifting polymers. 

1. Depinning of Interfaces 

1.1 Introduction and Phenomenology 

Depinning is a non-equilibrium critical phenomenon involving an external force 
and a pinning potential. When the force is weak the system is stationary, 
trapped in a metastable state. Beyond a threshold force the (last) metastable 
state disappears and the system starts to move. A simple example is provided 
by a point mass on a rough table. The mass is stationary until the external 
force F exceeds that of static friction F c . Larger forces lead to an initial period 
of acceleration, before the motion settles to a uniform velocity due to viscous 
forces. In the latter is proportional to velocity, the ultimate velocity of the 
point close to threshold behaves as v cx (F — F c ). 

While there are many other macroscopic mechanical examples, our main 
interest comes from condensed matter systems such as Charge Density Waves 
(CDWs) 1 , interfaces 2 , and contact lines 3 . In CDWs, the control parameter is 
the external voltage; a finite CDW current appears only beyond a threshold 
applied voltage. Interfaces in porous media, domain walls in random magnets, 
are stationary unless the applied force (magnetic field) is sufficiently strong. 
A key feature of these examples is that they involve the collective depinning of 
many degrees of freedom that are elastically coupled. As such, these problems 
belong to the realm of collective critical phenomena, characterized by universal 
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Figure 1. Geometry of the line in two dimensions. 



scaling laws. We shall introduce these laws and the corresponding exponents 
below for the depinning of a line (interface or contact line). 

Consider a line in two dimensions, oriented along the x direction, and 
fluctuating along a perpendicular direction r. The configuration of the line at 
time t is described by the function r(x,t). The function r is assumed to be 
single valued, thus excluding configurations with overhangs. In many cases 2 , 
where viscous forces dominate over inertia, the local velocity of a point on the 
curve is given by 

d ^=F + f{x,r)+K[r]. (1) 

The first term on the right hand side is a uniform applied force which is also the 
external control parameter. Fluctuations in the force due to randomness and 
impurities are represented by the second term. With the assumption that the 
medium is on average translationally invariant, the average of / can be set to 
zero. The final term in eq.(l)describes the elastic forces between different parts 
of the line. Short range interactions can be described by a gradient expansion; 
for example, a line tension leads to /C[r(a;)] = V 2 r or /C[r(g)] = —q 2 r(q) for 
the Fourier modes. The surface of a drop of non-wetting liquid terminates 
at a contact line on a solid substrate 3 . Deformations of the contact line are 
accompanied by distortions of the liquid/gas surface. As shown by Joanny 
and de Gennes 4 , the resulting energy and forces are non-local, described by 
JC[r(q)] = -\q\r(q). 

For the case of a surface in three dimensions deformations are described by 
r(xi,X2)- More generally, we shall consider r(x), where x is a d— dimensional 
vector. In a similar spirit, we shall generalize the coupling to /C[r((/)] = 
— \q\ a r(q), which interpolates between the above two cases as a changes from 
one to two. Note, however, that the equation of motion need not originate 
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from variations of a Hamiltonian, and may include non-linear couplings which 
will be discussed later on. 
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Figure 2. Critical behavior of the velocity. 



When F is small, the line is trapped in one of many metastable states in 
which dr/dt = at all points. For F larger than a threshold F c , the line is 
depinned from the last metastable state, and moves with an average velocity 
v. On approaching the threshold from above, the velocity vanishes as 



where j3 is the velocity exponent, and A is a nonuniversal amplitude. A mean- 
field estimate for (3 was obtained by Fisher in the context of CDWs 5 . It 
corresponds to the limit a = 0, where every point is coupled to all others, and 
hence experiences a restoring force proportional to (r(x)} — r(x). The resulting 
equation of motion, 



has to be supplemented with the condition (r(x)) = vt. The self-consistent 
solution for the velocity indeed vanishes as (F — Fc) 13 , with an exponent that 
depends on the details of the random force. If /(x, r(x)) varies smoothly with 
r, the exponent is /3 = 3/2, while discontinuous jumps in the force (like a 
saw-tooth) result in (3 = 1. In fact the latter is a better starting point for 
dcpinning in finite dimensions. This is because of the avalanches in motion 
(discussed next), which lead to a discontinuous coarse grained force. 

The motion just above threshold is not uniform, composed of rapid jumps 
as large segments of the line depin from strong pinning centers, superposed on 



v = A(F - F c f, 



(2) 



cir(x) 
dt 



= <r(x))-r(x)+F + /(x,r(x)), 
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the slower steady advance. These jumps have a power law distribution in size, 
cutoff at a correlation length £ which diverges at the transition as 

Z~(F-F C )- V . (3) 

The jumps are reminiscent of avalanches in other slowly driven systems. In 
fact, the depinning can be approached from below F c by monotonically in- 
creasing F in small increments, each sufficient to cause a jump to the next 
metastable state. The size and width of avalanches becomes invariant on ap- 
proaching F c . For example, 

Prob(width of avalanche > t) w ^p(£/£_), (4) 

where the cutoff diverges as in Eq.(3). The critical line is a self-affine 
fractal whose correlations satisfy the dynamic scaling from 

([r(x, t) r(x', t')f) = (x- X 'f C 9 ( j^Tjl) - ( 5 ) 

defining the roughness and dynamic exponents, £ and z respectively. (Angular 
brackets reflect averaging over all realizations of the random force /.) The 
scaling function g goes to a constant as its argument approaches 0; £ is the 
wandering exponent of an instantaneous line profile, and z relates the average 
lifetime of an avalanche to its size by r(£) <~ £ z . 

Although, the underlying issues of collective depinning for CDWs and in- 
terfaces have been around for some time, only recently a systematic perturba- 
tive approach to the problem was developed. This functional renormalization 
group (RG) approach to the dynamical equations of motion was originally de- 
veloped in the context of CDWs by Narayan and Fisher 6 (NF), and extended 
to interfaces by Nattermann et al 7 . We shall provide a brief outline of this ap- 
proach starting from Eq.(l). Before embarking on the details of the formalism, 
it is useful to point out some scaling relations amongst the exponents which 
follow from underlying symmetries and non-renormalization conditions. 

1. As mentioned earlier, the motion of the line close to the threshold is com- 
posed of jumps of segments of size £. Such jumps move the interface forward 
by ^ over a time period £ z . Thus the velocity behaves as, 

_jr c |"(*-0 = !/(*_£). (6) 

2. If the elastic couplings are linear, the response of the line to a static per- 
turbation e(x) is obtained simply by considering 

r £ (x,t)=r(x,t)-)C- 1 [e(x)}, (7) 
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where /C _1 is the inverse kernel. Since, r £ satisfies Eq.(l) subject to a force 
F + e(x) + f{x,r e ), r satisfies the same equation with a force F + f(x,r — 
JC~ 1 {e{x)\). As long as the statistical properties of the stochastic force are not 
modified by the above change in its argument, d (r)/de — 0, and 

ar <W\=-iC-\ or /5#\ = J_. (8) 



ds(x) I \ de(q) / \q 

Since it controls the macroscopic response of the line, the kernel JC cannot 
change under RG scaling. From Eqs.(5) and (3), we can read off the scaling 
of r(x), and the force SF, which using the above non-renormalization must be 
related by the exponent relation 

C+l=v. (9) 

Note that this identity depends on the statistical invariance of noise under 
the transformation in Eq.(7). It is satisfied as long as the force correlations 
(f(x,r)f(x',r')} only depend on r — r' . The identity does not hold if these 
correlations also depend on the slope dr/dx. 

3. A scaling argument related to the Imry-Ma estimate of the lower critical 
dimension of the random field Ising model, can be used to estimate the rough- 
ness exponent 8 . The elastic force on a segment of length £ scales as . If 
fluctuations in force are uncorrelated in space, they scale as £-« +1 )/ 2 over the 
area of an avalanche. Assuming that these two forces must be of the same 
order to initiate the avalanche leads to 

This last argument is not as rigorous as the previous two. Nonetheless, all 
three exponent identities can be established within the RG framework. Thus 
the only undetermined exponent is the dynamic one, z. 

1.2 Functional Renormalization Group 

A field theoretical description of the dynamics of Eq.(l) can be developed 
using the formalism of Martin, Siggia and Rose 9 (MSR): Generalizing to a 
d— dimensional interface, an auxiliary field f(x, t) is introduced to implement 
the equation of motion as a series of 5-functions. Various dynamical response 
and correlation functions for the field r(x, t) can then be generated from the 
functional, 

Z = J Vr(x,t)Vf(x,t)J[r]exp(S), (11) 
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where 

S = i J d d xdtf(x,t){d t r-JC[r]-F- f (x,r(x,t))}. (12) 

The Jacobian J[r] is introduced to ensure that the 5-functions integrate to 
unity. It does not generate any new relevant terms and will be ignored hence- 
forth. 

The disorder-averaged generating functional Z can be evaluated by a 
saddle-point expansion around a Mean-Field (MF) solution obtained by set- 
ting [r(x)] = vt — r(x). This amounts to replacing interaction forces with 
Hookean springs connected to the center of mass, which moves with a velocity 
v. The corresponding equation of motion is 

=vt- r MF {t) + f[r MF (t)] + F MF {v), (13) 

where the relationship Fmf{v) between the external force F and average ve- 
locity v is determined from the consistency condition (rMF(t)) = vt. The 
MF solution depends on the type of irregularity 6 : For smoothly varying ran- 
dom potentials, (3mf = 3/2, whereas for cusped random potentials, 0mf = 1- 
Following the treatment of NF 6 ' 10 , we use the mean field solution for cusped 
potentials, anticipating jumps with velocity of 0(1), in which case (3mf = 1- 
After rescaling and averaging over impurity configurations, we arrive at a gen- 
erating functional whose low-frequency form is 



Z = J VR(x,t)VR(x,t) exp(S), 

S = - J d d xdt[F- F MF (v)] R(x, t) 
r d d q duo 

~ J w 



(14) 

+ i J d d x dt dt' R(x, t)R(x, t')C [vt - vt' + R(x, i) - R(x, t')] . 

In the above expressions, R and R are coarse-grained forms of r — vt and if, re- 
spectively. F is adjusted to satisfy the condition (R) = 0. The function C(vt) 
is initially the connected mean-field correlation function {(rMF{t)rMF(t + T)) c . 

Ignoring the independent terms in the argument of C, the action becomes 
Gaussian, and is invariant under a scale transformation x — > bx, t — > 6 CT i, 
R -» b a - d / 2 R, R -» b-< y - d ' 2 R, F -> b- d ' 2 F, and v -» b- d / 2 v. Other terms 
in the action, of higher order in R and R, that result from the expansion of C 
[and other terms not explicitly shown in Eq.(14)], decay away at large length 
and time scales if d > d c = 2a. For d > d c , the interface is smooth (Co < 0) 



6 



at long length scales, and the depinning exponents take the Gaussian values 
z = o,vq = 2/d, 0o = 1- 

At d — d c , the action S has an infinite number of marginal terms that can 
be rearranged as a Taylor series for the function C [vt — vt' + i?(x, t) — i?(x, t')\, 
when v — > 0. The RG is carried out by integrating over a momentum shell 
A/6 < |q| < A (we set the cutoff wave vector to A = 1 for simplicity) and all 
frequencies, followed by a scale transformation x — > bx, t — ► 6 z t, i? — > and 
i? -» 6 e " d i?, where 6 = e*. The resulting recursion relation for the linear part 
in the effective action (to all orders in perturbation theory) is 

^ W M ^ ={Z + 9)i - F ~ Fmf) + constant ' (15) 
which immediately implies (with a suitable definition of F c ) 

d -^^=VF{F-F c ), (16) 
with the exponent identity 

y F = z + 6 = l/ V . (17) 

The functional renormalization of C(u) in d = 2a — e interface dimensions, 
computed to one-loop order, gives the recursion relation, 

^ e+2 , +2( _ )]cw + < „^ 

where Sd is the surface area of a unit sphere in d dimensions. NF showed 
that all higher order diagrams contribute to the renormalization of C as to- 
tal derivatives with respect to u, thus, integrating Eq.(18) at the fixed-point 
solution dC'*/d£ = 0, together with Eqs.(9) and (17), gives ( — e/3 to all 
orders in e, provided that J C* ^ 0. This gives Eq.(10) for a one-dimensional 
interface, as argued earlier. This is a consequence of the fact that C(u) re- 
mains short-ranged upon renormalization, implying the absence of anomalous 
contributions to (. 

The dynamical exponent z is calculated through the renormalization of p, 
the term proportional to Rd t R, which yields 

z = cr-2e/9 + 0(e 2 ), (19) 
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and using the exponent identity (6), 

(3 = 1 - 2e/9a + 0(e 2 ). (20) 

Nattermann et. al. obtain the same results to 0(e) by directly averaging the 
MSR generating function in Eq.(ll), and expanding perturbatively around a 
rigidly moving interface. 

Numerical integration of Eq.(l) for an elastic interface 11 (<r = 2) has 
yielded critical exponents £ = 0.97 ± 0.05 and v — 1.05 ± 0.1, in agreement 
with the theoretical result £ = v = 1. The velocity exponent (3 = 0.24 ± 0.1 is 
also consistent with the one-loop theoretical result 1/3; however, a logarithmic 
dependence v <~ 1/ ln(.F — F c ), which corresponds to j3 = 0, also describes the 
numerical data well. In contrast, experiments and various discrete models of 
interface growth have resulted in scaling behaviors that differ from system to 
system. A number of different experiments on fluid invasion in porous media 12 
give roughness exponents of around 0.8, while imbibition experiments 13 ' 14 have 
resulted in £ 0.6. A discrete model studied by Leschhorn 15 , motivated by 
Eq.(l) with cr = 2, gives a roughness exponent of 1.25 at threshold. Since the 
expansion leading to Eq.(l) breaks down when £ approaches one, it is not clear 
how to reconcile the results of Leschhhorn's numerical work 15 with the coarse- 
grained description of the RG calculation, especially since any model with 
£ > 1 cannot have a coarse grained description based on gradient expansions. 

1.3 Anisotropy 

Amaral, Barabasi, and Stanley (ABS) 16 recently pointed out that various mod- 
els of interface depinning in 1+1 dimensions fall into two distinct classes, de- 
pending on the tilt dependence of the interface velocity: 

1. For models like the random field Ising Model 17 , and some Solid On Solid 
models, the computed exponents are consistent with the exponents given by the 
RG analysis. It has been suggested 15 , however, that the roughness exponent 
is systematically larger than e/3, casting doubt on the exactness of the RG 
result. 

2. A number of different models, based on directed percolation (DP) 18 ' 13 give 
a different roughness exponent, £ « 0.63. In these models, pinning sites are 
randomly distributed with a probability p, which is linearly related to the force 
F. The interface is stopped by the boundary of a DP cluster of pinning sites. 
The critical exponents at depinning can then be related to the longitudinal 
and transverse correlation length exponents vu w 1.70 and v± w 1.07 of DP. 
In particular, £ = v\\/v±_ ~ 0.63, and (3 = v\\ — v± w 0.63, in agreement with 
experiments. 

The main difference of these models can be understood in terms of the 
dependence of the threshold force F c to the orientation. To include the possible 
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dependence of the line mobility on its slope, d x r, we can generalize the equation 
of motion to 

t r = Kd 2 x r + Kd x r + ^(d x rf + F + f(x, r). (21) 

The isotropic depinning studied by RG corresponds to k = A = 0. The usual 
mechanisms for generating a non-zero A are of kinematic origin 19 (A oc v) and 
can be shown to be irrelevant at the depinning threshold where the velocity 
v goes to zero 10 . However, if A is not proportional to v and stays finite at 
the transition, it is a relevant operator and expected to modify the critical 
behavior. As we shall argue below, anisotropy in the medium is a possible 
source of the nonlinearity at the depinning transition. 

A model flux line (FL) confined to move in a plane 11 ' 20 provides an ex- 
ample where both mechanisms for the nonlinearity are present. Only the force 
normal to the FL is responsible for motion, and is composed of three com- 
ponents: (1) A term proportional to curvature arising from the smoothening 
effects of line tension. (2) The Lorentz force due to a uniform current density 
perpendicular to the plane acts in the normal direction and has a uniform mag- 
nitude F (per unit line length). (3) A random force n • f due to impurities, 
where n is the unit normal vector 20 . Equating viscous dissipation with the 
work done by the normal force leads to the equation of motion 

where h(x,t) denotes transverse displacement of the line and s = d x h. The 
nonlinearities generated by \f\ + s 2 are kinematic in origin 19 and irrelevant as 
v — > 10 , as can be seen easily by taking them to the left hand side of Eq.(22). 
The shape of the pinned FL is determined by the competition of the terms in 
the square brackets. Although there is no explicit simple s 2 term in this group, 
it will be generated if the system is anisotropic. 

To illustrate the idea, let us take fh and f x to be independent random 

1/2 1/2 

fields with amplitudes A ft and A x respectively; each correlated isotropically 
in space within a distance a. For weak disorder, a deformation of order a in 
the normal direction n takes place over a distance L c 3> a along the line. The 
total force due to curvature on this piece of the line is of the order of L c (a/L 2 ), 
and the pinning force, [(L c /a)(n 2 A/, + n 2 A^)] 1 / 2 . Equating the two forces 2 
yields L c = a(n\Ah + n x A^) -1 / 3 and an effective pinning strength per unit 
length, 

9 1 /A h + s 2 A a; \2/3 
F ( 8 ) = oL e - 2 =o- 1 ( \ + g2 x ) . 

The roughening by impurities thus reduces the effective driving force on the 
scale L c to F(s) = F — Fq(s). Therefore, even if initially F is independent 



dih 



(1 + s 2 )3/2 



fh ~ Sfx 

VTT^ 2 



(22) 
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of s, such a dependence is generated under coarse graining, provided that the 
random force is anisotropic, i.e. 7^ A^. An expansion of F(s) around its 
maximum (which defines the hard direction) yields an s 2 term which is positive 
and remains finite as v — > 0. 

The above example indicates the origin of the two types of behavior for 
A ff = v"(s = 0) observed by ABS 16 : Kinematics produces a A e ff proportional 
to v which vanishes at the threshold; anisotropy yields a nonvanishing (and 
diverging) A e fj at the depinning transition. An immediate consequence of the 
latter is that the depinning threshold F c depends on the average orientation 
of the line. While anisotropy may generate other local terms in the effective 
equation of motion, at a symmetry direction, this term is the only relevant one 
in the RG sense, capable of modifying the critical behavior for d < 4. A one- 
loop RG of Eq.(21) with the k = was carried out by Stepanow 21 . He finds 
no stable fixed point for 2 < d < 4, but his numerical integration of the one 
loop RG equations in d = 1 yield ( 0.8615 and a dynamical exponent z = \. 
Due to the absence of Galilean invariance, there is also a renormalization of A 
which is related to the diverging A e gf observed in Ref. 16 . The nonperturbative 
nature of the fixed point precludes a gauge of the reliability of these exponents. 

Numerical simulations of Eq.(21) in d = l 22 , indicate that it shares the 
characteristics of a class of lattice models 18 ' 13 where the external force is re- 
lated to the density p of "blocking sites" by F = 1 — p. When p exceeds a 
critical value of p c , blocking sites form a directed percolating path which stops 
the interface. For a given geometry, there is a direction along which the first 
spanning path appears. This defines a hard direction for depinning where the 
threshold force F c (s) reaches maximum. Higher densities of blocking sites are 
needed to form a spanning path away from this direction, resulting in a lower 
threshold force F c (s) for a tilted interface. Thus on a phenomenological level 
we believe that the nonlinear equation, and directed percolation (DP) models 
of interface depinning belong to the same universality class of anisotropic de- 
pinning. This analogy may in fact be generalized to higher dimensions, where 
the blocking path is replaced by a directed blocking surface 23 . Unfortunately, 
little is known analytically about the scaling properties of such a surface at 
the percolation threshold. 

As emphasized above, the hallmark of anisotropic depinning is the de- 
pendence of the threshold force F c (s) on the slope s. Above this threshold, 
we expect v(F, s) to be an analytical function of F and s. In particular, for 
F > F c (0), there is a small s expansion v(F, s) — v(F, s = 0)+A o ffS 2 /2+- • •. On 
the other hand, we can associate a characteristic slope s = ~ (SF)"^- 1- ^ , 

to DP clusters where SF — F — F c (0), and v is the correlation length exponent. 
Scaling then suggests 

v(F,s) = (8F) 6 g(s/SF^ 1 -^), (23) 
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where 6 = v(z — (). Matching Eq.(23) with the small s expansion, we see 
that Aoff diverges as (dF)-* (as defined by ABS 16 ) with = 2v{\ - Q - 6 = 
u(2 — ( — z). In d — 1, the exponents v and £ are related to the correlation 
length exponents v\\ and of DP 23 via v — vn k, 1.73 and ( = v_\_/v\\ w 0.63, 
while the dynamical exponent is z = 1. Scaling thus predicts 4> w 0.63, in 
agreement with the numerical result of 0.64 ± 0.08 in Ref. 16 . Close to the line 
F = F c (0) (but at a finite s), the dependence of v on SF drops out and we 
have 

v(F c ,s) cx \s\ e ^ (1 - c \ (24) 

As2 = lind=l, the above equation reduces to v cx \s\, in agreement with 
Fig. 1 of Ref. 16 . Since v(F, s) = at F = F c (s), Eq.(23) suggests 

FcW-Fc^a-l*! 1 /"^. (25) 

Note that Eqs. (24) and (25) are valid also in higher dimensions, though values 
of the exponents quoted above vary with d 23 . 

An interface tilted away from the hard direction not only has a different 
depinning threshold, but also completely different scaling behavior at its tran- 
sition. This is because, due to the presence of an average interface gradient 
s = (V/i), the isotropy in the internal x space is lost. The equation of motion 
for fluctuations, /i'(x, t) — /i(x, t) — s • x, around the average interface position 
may thus include a non-zero k in (21). The resulting depinning transition be- 
longs to yet a new universality class with anisotropic response and correlation 
functions in directions parallel and perpendicular to s; i.e. 



([Mx)-M^)] 2 )=k 




where n is the ansiotropy exponent, and x t denotes the d—1 directions trans- 
verse to s. 

A suggestive mapping allows us to determine the exponents for depinning 
a tilted interface: Consider the response to a perturbation in which all points 
along a (d— l)-dimcnsional cross section of the interface at a fixed x\\ are pushed 
up by a small amount. This move decreases the slope of the interface uphill 
but increases it downhill. Since F c (s) decreases with increasing s, at criticality 
the perturbation propagates only a finite distance uphill but causes a downhill 
avalanche. The disturbance front moves at a constant velocity (5xu cx t) and 
hence 211 =1. (Such chains of moving sites were indeed seen in simulations of 
the d = 2 model discussed below.) Furthermore, the evolution of successive 
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cross sections x t (x||) is expected to be the same as the evolution in time of 
a (d — l)-dimensional interface! The latter is governed by the Kardar-Parisi- 
Zhang (KPZ) equation 19 , whose scaling behavior has been extensively studied. 
From this analogy we conclude, 

ZKPz{d-l) z K p Z {d-l) 

In particular, the tilted interface with d = 2 maps to the growth problem in 
1+1 dimensions where the exponents are known exactly, yielding £(2) = 1/3 
and 17(2) = 2/3. This picture can be made more precise for a lattice model 
introduced below. Details will be presented elsewhere. 

To get the exponent (3 for the vanishing of velocity of the tilted interface, 
we note that since zu = 1, v scales as the excess slope 5s = s — s c (F). The 
latter controls the density of the above moving fronts; s c (F) is the slope of 
the critical interface at a given driving force F, i.e., F — F c (s c ). Away from 
the symmetry direction, the function F c (s) has a non- vanishing derivative and 
hence 

6F = F- F c (s) = F c {s c ) - F c {s) ~ 6s ~ v. (27) 

We thus conclude that generically (3 = 1 for tilted interfaces, independent of 
dimension. 

To check the above predictions, we performed simulations of the paral- 
lelized version of a previously studied percolation model of interface depinning 18 . 
A solid-on-solid (SOS) interface is described by a set of integer heights {hi} 
where i is a group of d integers. With each configuration is associated a ran- 
dom set of pinning forces {iji G [0,1)}. The heights are updated in parallel 
according to the following rules: hi is increased by one if (i) hi < h- } — 2 for 
at least one j which is a nearest neighbor of i, or (ii) rji < F for a pre-selected 
uniform force F. If hi is increased, the associated random force rji is also up- 
dated, i.e. replaced by a new random number in the interval [0, 1). Otherwise, 
hi and rji are unchanged. The simulation is started with initial conditions 
hi(t = 0) = Int[siJ, and boundary conditions /m+l = Int [sL] + hi are enforced 
throughout. The CPU time is greatly reduced by only keeping track of active 
sites. 

The above model has a simple analogy to a resistor-diode percolation 
problem 23 . Condition (i) ensures that, once a site (i, h) is wet (i.e., on or 
behind the interface), all neighboring columns of i must be wet up to height 
h — 1. Thus there is always "conduction" from a site at height h to sites in 
the neighboring columns at height h — 1. This relation can be represented 
by diodes pointing diagonally downward. Condition (ii) implies that "conduc- 
tion" may also occur upward. Hence a fraction F of vertical bonds are turned 
into resistors which allow for two-way conduction. Note that, due to the SOS 
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condition, vertical downward conduction is always possible. For F < F c , con- 
ducting sites connected to a point lead at the origin, form a cone whose hull 
is the interface separating wet and dry regions. The opening angle of the cone 
increases with F, reaching 180° at F = F c , beyond which percolation in the 
entire space takes place, so that all sites are eventually wet. If instead of a 
point, we start with a planar lead defining the initial surface, the percolation 
threshold depends on the surface orientation, with the highest threshold for 
the untiltcd one. 

Our simulations of lattices of 65536 sites in d = 1 and of 512 x 512 and 
840 x 840 sites in d = 2 confirm the exponents for depinning in the hard 
direction. For a tilted surface in d = 1 the roughness exponent determined from 
the height-height correlation function is consistent with the predicted value of 
£ = 1/2 and different from ( w 0.63 of the untiltcd one. The dependence 
of the depinning threshold on slope is clearly seen in the figure below, where 
the average velocity is plotted against the driving force for s = (open) and 
s = 1/2 (solid). The s = data can be fitted to a power-law v ~ (F — F c ) 6 , 
where F c w 0.461, (3 = 0.63 ± 0.04 for d = 1, and F c w 0.201, (3 = 0.72 ± 0.04 
for d—2. Data at s = 1/2 are consistent with Eq.(27) close to the threshold. 

.35 
.30 
.35 
.80 
. 15 
. 10 
.05 


.1 .2 .3 .4 .5 .6 

F 

Figure 3. Average interface velocity v versus the driving force F, 
for d = 1, s = (open circles), d = 1, s = 1/2 (solid circles), d = 2, 
s = (open squares), and d~ 2, s = 1/2 (solid squares). 

We also measured height-height correlation functions at the depinning 
transition. For a tilted surface in d = 2, the height fluctuations and correspond- 
ing dynamic behaviors are different parallel and transverse to the tilt. The next 
figure shows a scaling plot of (a) C\\ (r« , t) = ([/ifei + rii , Xt, t) — h(xii , x t , t)] 2 ) 
and (b) C t {r t ,t) = ([h(x\\ , x t + r t , t) — h{x\\ , x t , t)] 2 } against the scaled distances 
at the depinning threshold of an s = 1/2 interface. Each curve shows data at 
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a given t = 32, 64, • • •, 1024, averaged over 50 realizations of the disorder. The 
data collapse is in agreement with the mapping to the KPZ equation in one 
less dimension. 




r„ A r t /t z/3 

Figure 4. Height-height correlation functions (a) along and (b) 
transverse to the tilt for an 840 2 system at different times 32 < t < 
1024. The interface at t = is flat; d = 2, s = 1/2, and F = 0.144. 



In summary, critical behavior at the depinning of an interface depends on 
the symmetries of the underlying medium. Different universality classes can be 
distinguished from the dependence of the threshold force (or velocity) on the 
slope, which is reminiscent of similar dependence in a model of resistor-diode 
percolation. In addition to isotropic depinning, we have so far identified two 
classes of anisotropic depinning: along a (hard) axis of inversion symmetry in 
the plane, and tilted away from it. We have no analytical results in the former 
case, but suggest a number of scaling relations that are validated by simula- 
tions. In the latter (more generic) case we have obtained exact information 
from a mapping to moving interfaces, and confirmed them by simulations in 
d = 1 and d = 2. As it is quite common to encounter (intrinsic or artificially 
fabricated) anisotropy for flux lines in superconductors, domain walls in mag- 
nets, and interfaces in porous media, we expect our results to have important 
experimental ramifications. 

Another form of anisotropy is also possible for interfaces in 2+1 dimen- 
sions. If the directions x and y on the surface are not related by symmetry, 
the non-linear term in the KPZ equation can be generalized, leading to the 
depinning equation 

d t h = K x d 2 x r + K y dlr + ^(d x rf + ^(d y rf + F + f(x,y, r). (28) 
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In fact the difference between K x and K y is not important as long as both 
are positive. It was first pointed out by Dietrich Wolf 24 that different signs 
of X x and X y lead to a different universality class for the case of annealed 
noise. More recently it was demonstrated by Jeong et al 25 that, with quenched 
noise, eq.(28)describes a new universality class of depinning transitions with 
w 0.80(1), and anisotropic roughness exponents in the x and y directions. 

2 Fluctuating Lines 

2.1 Flux Line Depinning 

The pinning of flux lines (FLs) in Type-II superconductors is of fundamen- 
tal importance to many technological applications that require large critical 
currents 26 . Upon application of an external current density J, the motion of 
FLs due to the Lorentz force causes undesirable dissipation of supercurrents. 
Major increases in the critical current density J c of a sample are achieved 
when the FLs are pinned to impurities. There are many recent studies, both 
experimental 27,28 and theoretical 29,30 , on collective pinning of FL's to point or 
columnar defects. Another consequence of impurities is the strongly nonlinear 
behavior of the current slightly above the depinning threshold, as the FLs start 
to move across the sample. Recent numerical simulations have concentrated 
on the low temperature behavior of a single FL near depinning 31,11,20 , mostly 
ignoring fluctuations transverse to the plane defined by the magnetic field and 
the Lorentz force. Common signatures of the depinning transition from J < J c 
to J > J c include a broad band (f~ a type) voltage noise spectrum, and self- 
similar fluctuations of the FL profile. 

The FL provides yet another example of a depinning transition. We now 
extend the methods of the previous section to the full three-dimensional dy- 
namics of a single FL at low temperatures. The shape of the FL at a given 
time t is described by r(x,t), where x is along the magnetic field B, and the 
unit vector is along the Lorentz force F. Point impurities are modeled by 
a random potential V(x, r), with zero mean and short-range correlations. In 
the presence of impurities and a bulk Lorentz force F, the energy of a FL with 
small fluctuations is, 



The simplest possible Langevin equation for the FL, consistent with local, 
dissipative dynamics, is 



H = 




(29) 



(30) 
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Figure 5. Geometry of the line in three dimensions. 



where fi is the mobility of the FL, and f = — V r V. The potential V(x,r) need 
not be isotropic. For example, in a single crystal of ceramic superconductors 
with the field along the oxide planes, it will be easier to move the FL along 
the planes. This leads to a pinning threshold that depends on the orientation 
of the force. Anisotropy also modifies the line tension, and the elastic term in 
Eq.(30) is in general multiplied by a non-diagonal matrix K a p. The random 
force f (x, r), can be taken to have zero mean with correlations 



We shall focus mostly on the isotropic case, with A Q7 (r — r') = 5 a7 A(|r — r'|), 
where A is a function that decays rapidly for large values of its argument. 

While the flux line is pinned by impurities when F < F c , for F slightly 
above threshold, we expect the average velocity v — |v| to scale as in Eq.(23). 
Superposed on the steady advance of the FL are rapid "jumps" as portions of 
the line depin from strong pinning centers. The cut off length £ on avalanche 
sizes diverges on approaching the threshold as £ <~ (F — F c )~ v . At length scales 
up to £, the correlated fluctuations satisfy the dynamic scaling forms, 



where £ Q and z a are the roughness and dynamic exponents, respectively. The 
scaling functions g a go to a constant as their arguments approach 0. Beyond 



(f a (x,r)f 1 (x',r' 



))=S(x-x')A aj (v-r'). 



(31) 



<[r||(s,t)-r||(0,0)] 2 ) =\x\ 2 ^ gil (t/\xn, 
(K(:M)-rx(0,0)] 2 ) =\ x \^ g± (t/\xn 



(32) 
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the length scale £, different regions of the FL dcpin more or less independently 
and the system crosses over to a moving state, described by different exponents, 
which will be considered in the next section. 

The major difference of this model from the previously studied interface is 
that the position of the flux line, r(x, t), is now a 2-dimensional vector instead 
of a scalar; fluctuating along both e\\ and directions. One consequence is 
that a "no passing" rule 32 , applicable to CDWs and interfaces, does not ap- 
ply to FLs. ft is possible to have coexistence of moving and stationary FLs 
in particular realizations of the random potential. How do these transverse 
fluctuations scale near the depinning transition, and do they in turn influ- 
ence the critical dynamics of longitudinal fluctuations near threshold? The 
answer to the second question can be obtained by the following qualitative 
argument: Consider Eq.(30) for a particular realization of randomness f(x,r). 
Assuming that portions of the FL always move in the forward direction, there 
is a unique point r±_{x,r\\) that is visited by the line for given coordinates 
{x,r\\). We construct a new force field /' on a two dimensional space {x,r\\) 
through f'(x, r|| ) = /|| (x,r||,r^(x,r||)). ft is then clear that the dynamics of 
the longitudinal component r\\{x,t) in a given force field f(x,r) is identical 
to the dynamics of r\\{x,t) in a force field f'(x, 7*||), with r± set to zero. It 
is quite plausible that, after averaging over all f, the correlations in /' will 
also be short-ranged, albeit different from those of f. Thus, the scaling of 
longitudinal fluctuations of the depinning FL will not change upon taking into 
account transverse fluctuations. However, the question of how these transverse 
fluctuations scale still remains. 

Certain statistical symmetries of the system restrict the form of response 
and correlation functions. For example, Eq.(30) has statistical space- and time- 
translational invariance, which enables us to work in Fourier space, i.e. (x, i) — > 
(q, lo). For an isotropic medium, F and v are parallel to each other, i.e., v(F) = 
v(F)F, where F is the unit vector along F. Furthermore, all expectation values 
involving odd powers of a transverse component are identically zero due to 
the statistical invariance under the transformation r± — ► — r±. Thus, linear 
response and two-point correlation functions are diagonal. The introduced 
critical exponents are then related through scaling identities. These can be 
derived from the linear response to an infinitesimal external force field e(q, u), 

Xap(q,u) = = S aPXa{q,u), (33) 

in the (q,oj) — > (0,0) limit. Eq.(30) is statistically invariant under the trans- 
formation F — ► F + e(q), r(q,w) — > r(q,w) + q~ 2 e(q). Thus, the static linear 
response has the form X||(s> w = 0) = xj_(g,a; = 0) = q~ 2 . Since £|| scales like 
the applied force, the form of the linear response at the correlation length £ 
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gives the exponent identity 

C|| + 1/f = 2. (34) 

Considering the transverse linear response seems to imply (± — Q . How- 
ever, the static part of the transverse linear response is irrelevant at the critical 
RG fixed point, since zj_ > z\\ , as shown below. When a slowly varying uniform 
external force e(t) is applied, the FL responds as if the instantaneous external 
force F + e is a constant, acquiring an average velocity, 



(d t r a ) = v a {F + e)^ v a (F) + ^r£ 7 - 



Substituting dv\\/dF\\ = dv/dF and dv±/dF± = v/F, and Fourier transform- 
ing, gives 

X\\(q = 0,w) = . — — j 

-I (35) 

X±(q = 0,w) 



' -ilo{v / F)- 1 + 0{lo 2 Y 

Combining these with the static response, we see that the characteristic relax- 
ation times of fluctuations with wavelength £ are 



dv^ 1 



r±(« = r 1 )~(fl a J)" 1 ~^ /,/, '~€* i , 

which, using Eq.(34), yield the scaling relations 

/? = (*u-CiiK 

Zj_ = Zn + \ jv. 



(36) 



We already see that the dynamic relaxation of transverse fluctuations is much 
slower than longitudinal ones. All critical exponents can be calculated from 
Oi, C-Lj an d by using Eqs(34), and (36). 

Equation (30) can be analyzed using the formalism of Martin, Siggia, 
and Rose (MSR) 9 . Ignoring transverse fluctuations, and generalizing to d 
dimensional internal coordinates xe5f d , leads to an interface depinning model 
which was studied by Nattermann, Stepanow, Tang, and Leschhorn (NSTL) 7 , 
and by Narayan and Fisher (NF) 10 . The RG treatment indicates that impurity 
disorder becomes relevant for dimensions d < 4, and the critical exponents in 
d = 4— e dimensions are given to one-loop order as ( — e/3, z = 2 — 2e/9. NSTL 
obtained this result by directly averaging the MSR generating functional Z, 
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and calculating the renormalization of the force-force correlation function A(r), 
perturbatively around the freely moving interface [A(r) = 0]. NF, on the other 
hand, used a perturbative expansion of around a saddle point corresponding 
to a mean- field approximation to Eq.(30) 33 , which involved temporal force- force 
correlations C(vt). They argue that a conventional low- frequency analysis 
is not sufficient to determine critical exponents. They also suggest that the 
roughness exponent is equal to e/3 to all orders in perturbation theory. 

Following the approach of NF, we employ a perturbative expansion of 
the disorder-averaged MSR partition function around a mean-field solution for 
cusped impurity potentials 10 . All terms in the expansion involving longitudinal 
fluctuations are identical to the interface case, thus we obtain the same critical 
exponents for longitudinal fluctuations, i.e., £|| = e/3, z\\ = 2 — 2e/9 + 0(e 2 ). 
Furthermore, for isotropic potentials, the renormalization of transverse tem- 
poral force- force correlations C±(vt) yields a transverse roughness exponent 
(± = 5(||/2 — 2, to all orders in perturbation theory. For the FL (e = 3), the 
critical exponents are then given by 

C||=l, Z||«4/3, i/=l, 
13 « 1/3, C± - 1/2, z± « 7/3. 




Figure 6. A plot of average velocity versus external force for a sys- 
tem of 2048 points. Statistical errors are smaller than symbol sizes. 
Both fits have three adjustable parameters: The threshold force, the 
exponent, and an overall multiplicative constant. 



To test the scaling forms and exponents predicted by Eqs.(23) and (32), 
we numerically integrated Eq.(30), discretized in coordinates x and t. Free 
boundary conditions were used for system sizes of up to 2048, with a grid 
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Ix-x'l 

Figure 7. A plot of equal time correlation functions versus sepa- 
ration, for the system shown in Fig. 6, at F = 0.95. The observed 
roughness exponents very closely follow the theoretical predictions 
of £|| = 1, C-L = 0-5, which are shown as solid lines for comparison. 



spacing Ax = 1 and a time step At = 0.02. Time averages were evaluated 
after the system reached steady state. Periodic boundary conditions gave sim- 
ilar results, but with larger finite size effects. Smaller grid sizes did not change 
the results considerably. The behavior of v(F) seems to fit the scaling form 
of Eq.(23) with an exponent (3 s=s 0.3, but is also consistent with a logarithmic 
dependence on the reduced force, i.e., (i = 0. The same behavior was observed 
by Dong et al. in a recent simulation of the 1 + 1 dimensional geometry 11 . Since 
Z|| , and consequently f3, is known only to first order in e, higher order correc- 
tions are expected. By looking at equal time correlation functions, we find 
that transverse fluctuations are strongly suppressed, and that the roughness 
exponents are equal to our theoretical estimates within statistical accuracy. 
The excellent agreement for e = 3 suggests that the theoretical estimates are 
indeed exact. 

The potential pinning the FL in a single superconducting crystal is likely 
to be highly anisotropic. For example, consider a magnetic field parallel to 
the copper oxide planes of a ceramic superconductor. The threshold force 
then depends on its orientation, with depinning easiest along the copper oxide 
planes. In general, the average velocity may depend on the orientations of 
the external force and the FL. The most general gradient expansion for the 
equation of motion is then, 

^ = VapFfj + K a f3d x r fj + K aj3 dlr p + \k cttj3l d x r l3 d x r 1 + f a (x, r(x, t)) H , 
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with 



(f a (x, v)f p (x',r')) = S(x - x')C aP {v - r'). 



(39) 



Depending on the presence or absence of various terms allowed by the sym- 
metries of the system, the above set of equations encompasses many distinct 
universality classes. For example, consider the situation where v depends on 
F, but not on the orientation of the line. Eqs.(35) have to be modified, since v 
and F arc no longer parallel (except along the axes with r — ► — r symmetry), 
and the linear response function is not diagonal. The RG analysis is more 
cumbersome: For depinning along a non-symmetric direction, the longitudinal 
exponents are not modified (in agreement with the argument presented ear- 
lier), while the transverse fluctuations are further suppressed to £j_ = 2G| — 2 
(equal to zero for £y = l) 34 . Relaxation of transverse modes are still character- 
ized by z± — Z)| + 1/v, and the exponent identity (34) also holds. Surprisingly, 
the exponents for depinning along axes of reflection symmetry are the same as 
the isotropic case. 

If the velocity also depends on the tilt, there will be additional relevant 
terms in the MSR partition function, which invalidate the arguments leading 
to Eq.(34). The analogy to FLs in a planes suggests that the longitudinal 
exponents for d = 1 are controlled by DP clusters 18,13 , with Q w 0.63. Since 
no perturbativc fixed point is present in this case, it is not clear how to explore 
the behavior of transverse fluctuations systematically. 

2.2 Dynamic Fluctuations of an Unpinned Flux Line 

So far,w investigated the dynamics of a Flux Line near the depinning 
transition. Now, we would like to consider its behavior in a different regime, 
when the external driving force is large, and the impurities appear as weak 
barriers that deflect portions of the line without impeding its overall drift. 
In such non-equilibrium systems, one can regard the evolution equations as 
more fundamental, and proceed by constructing the most general equations 
consistent with the symmetries and conservation laws of the situation under 
study 35 . Even in a system with isotropic randomness, which we will discuss 
here, the average drift velocity, v, breaks the symmetry between forward and 
backward motions, and allows introduction of nonlinearities in the equations 
of motion 36 ' 35 . 

Let us first concentrate on an interface in two dimensions. (Fig.l.) By 
contracting up to two spatial derivatives of r, and keeping terms that are 
relevant, one obtains the Kardar-Parisi-Zhang 19 (KPZ) equation, 



d t r(x, t)=nF + Kd 2 x r{x, t) + - [d x r(x, t)f + f(x, t) 



(40) 



with random force correlations 



{f{x, t)f(x', t 1 )) = 2TS(x - x')5(t - t'). 



(41) 
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For a moving line, the term proportional to the external force can be absorbed 
without loss of generality by considering a suitable Galilean transformation, 
r — > r — at, to a moving frame. A large number of stochastic nonequilibrium 
growth models, like the Eden Model and various ballistic deposition models are 
known to be well described, at large length scales and times, by this equation, 
which is intimately related to several other problems. For example, the trans- 
formation v(x, t) = —\d x r(x, t) maps Eq.(40) to the randomly stirred Burgers' 
equation for fluid flow 37,38 , 

d t v + vd x v = Kd 2 x v - Xd x f{x, t). (42) 

The correlations of the line profile still satisfy the dynamic scaling form in 
Eq.(5), nevertheless with different scaling exponents £, z and scaling function 
g. This sclf-affine scaling is not critical, i.e., not obtained by fine tuning an 
external parameter like the force, and is quite different in nature than the 
critical scaling of the line near the depinning transition, which ceases beyond 
the correlation length scale £. 

Two important nonperturbative properties of Eq.(40) help us determine 
these exponents exactly in 1+1 dimensions: 

1. Galilean Invariance (GI): Eq.(40) is statistically invariant under the in- 
finitesimal reparametrization 

r' = r + ex , x 1 = x + Xet , t' = t, (43) 

provided that the random force / does not have temporal correlations 39 . Since 
the parameter A appears both in the transformation and Eq.(40), it is not 
rcnormalizcd under any RG procedure that preserves this invariance. This 
implies the exponent identity 38 ' 39 

C + z = 2. (44) 

2. Fluctuation-Dissipation (FD) Theorem: Eqs.(40) and (41) lead to aFokkcr 
Planck equation for the evolution of the joint probability V [r(x)] , 

f ( 5V S 2 V \ 

d ^ = j dx {m dtr + T ¥w)- (45) 

It is easy to check that V has a stationary solution 

V = cxp(-^ J dx(d x r) 2 y (46) 

If V converges to this solution, the long-time behavior of the correlation func- 
tions in Eq.(5) can be directly read off Eq.(46), giving ( = 1/2. 
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Combining these two results, the roughness and dynamic exponents are 
exactly determined for the line in two dimensions as 



C = 1/2 , z = 3/2. (47) 

Many direct numerical simulations and discrete growth models have verified 
these exponents to a very good accuracy. Exact exponents for the isotropic 
KPZ equation are not known in higher dimensions, since the FD property is 
only valid in two dimensions. These results have been summarized in a number 
of recent reviews 40 ' 41 ' 42 ' 43 . 

As an aside we remark that some exact information is available for the 
anisotropic KPZ equation in 2+1 dimensions. Using a perturbative RG ap- 
proach, Wolf showed 24 that in the equation 

d t r = KV 2 r + ^{d x rf + ^(d y rf + f(x, y, t), (48) 

the nonlinearities {X x ,X y } renormalize to zero if they initially have opposite 
signs. This suggests logarithmic fluctuations for the resulting interface, as 
in the case of the linear Langevin equation. In fact, it is straightforward to 
demonstrate that eq.(48) also satisfies a Fluctuation Dissipation condition if 
X x = —X y . When this condition is satisfied, the associated Fokker-Planck 
equation has a steady state solution 



(49) 



This is a non-perturbative result which again indicates the logarithmic fluc- 
tuations resulting from eq.(48)In this context, it is interesting to note that 
the steady state distribution for an exactly solvable discrete model of surface 
growth belonging to the above universality class has also been obtained 44 

Let us now turn to the case of a line in three dimensions (Fig. 5). Fluctu- 
ations of the line can be indicated by a a two dimensional vector r. Even in an 
isotropic medium, the drift velocity v breaks the isotropy in r by selecting a 
direction. A gradient expansion up to second order for the equation of motion 
gives 45 

d t r a = [Ki5 a p + K 2 v a v p ] d^rp 

d T rQ t ' ^^^^ 

+ [Ai(£ a/3 u 7 + 5 ai Vf3) + X 2 v a 8 Pl + X 3 v a vpv 7 ] - x 7 + f a 

with random force correlations 

(f a (x,t)ff 3 (x',t')) = 2[T 1 6 af3 + T 2 v a Vf 3 }6(x-x , )5(t-t'). (51) 
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Higher order nonlinearities can be similarly constructed but are in fact irrele- 
vant. In terms of components parallel and perpendicular to the velocity, the 
equations are 



d n = + ^{d xn f + ^(d x r ± f + /n(*,t) 

d t r ± = K±dlr± + X±d x r»d x r± + f±(x,t) 



(52) 



with 



{f\\{x,t)f\\{x',t')) =2T\\5{x-x')5{t-t') 
(f ± (x,t)f ± (x',t')) =2T ± 6(x - x')S(t - t') 
The noise-averaged correlations have a dynamic scaling form like Eq.(32) 

([ n {x,t)~ n (x',t>)?) = \x-x'\^g/ " ' 
([r ± (x,t)-r ± (x',t')} 2 ) = \x-xf^g ± 



(53) 



- x'\ z ± 



\X — X' 

In the absence of nonlinearities (A|| = A x = Ax = 0), Eqs.(52) can easily 
be solved to give Cj| — (± — 1/2 and zu = z± = 2. Simple dimensional 
counting indicates that all three nonlinear terms are relevant and may modify 
the exponents in Eq.(54). Studies of related stochastic equations 46,24 indicate 
that interesting dynamic phase diagrams may emerge from the competition 
between nonlinearities. Let us assume that A|| is positive and finite (its sign 
can be changed by rii — > — rn), and focus on the dependence of the scaling 
exponents on the ratios Aj_/A|| and A x /A||, as depicted in Fig. 8. (It is more 
convenient to set the vertical axis to A x -?^||Tj_/A||-ff_i_T||.) 

The properties discussed for the KPZ equation can be extended to this 
higher dimensional case: 

1. Galilean Invariance (GI): Consider the infinitesimal rcparamctrization 

x + \i\et , t' = t , 

11 , (55) 

= rn + ex , r± = rj_ . 



Eqs.(52) are invariant under this transformation provided that A|| = Aj_. Thus 
along this line in Fig. 8 there is GI, which once more implies the exponent 
identity 

C||+*||=2. (56) 
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A„K„T, 




Figure 8. A projection of RG flows in the parameter space, for n = 1 
transverse components. 



2. Fluctuation-Dissipation (FD) Condition: The Fokker-Planck equation for 
the evolution of the joint probability V [r\\ (x) , r± (x)] has a stationary solution 



P oc exp 



(-/ 



dx 



§(^n) a + ^(W 



(57) 



provided that A x K»T± = X±K±T^. Thus for this special choice of parameters, 
depicted by a starred line in Fig. 8, if V converges to this solution, the long- 
time behavior of the correlation functions in Eq.(54) can be directly read off 
Eq.(57), giving C|| = C± = 1/2- 

3. The Cole-Hopf ( CH) Transformation is an important method for the exact 
study of solutions of the one component nonlinear diffusion equation 37 . Here 
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we generalize this transformation to the complex plane by defining, for A x < 0, 



,w ,x ( \n{m,t) + iy/-\\\X x r±(x,t)\ 

*(x, t) = exp |- . (58) 



The linear diffusion equation 



then leads to Eqs.(52) if ^ = K ± = K and A|| = A ± . [Here Re(/i) = A||/||/2iT 
and Im(/z) = y 7 — A~jj A x f±/2K.] This transformation enables an exact solution 
of the deterministic equation, and further allows us to write the solution to 
the stochastic equation in the form of a path integral 



p(x,t) f ft 

= / Vx{t) exp < — / dr 
J(o,o) I 



x 2 



(59) 



Eq. (59) has been extensively studied in connection with quantum tunneling in 
a disordered medium 47 , with ^ representing the wave function. In particular, 
results for the tunneling probability |^| 2 suggest z\\ =3/2 and C\\ = 1/2. The 
transverse fluctuations correspond to the phase in the quantum problem which 
is not an observable. Hence this mapping does not provide any information on 
(± and z± which are in fact observable for the moving line. 

At the point A^ = A x = 0, ru and r± decouple, and z± = 2 while zu = 3/2. 
However, in general z\\ = z± = z unless the effective A^ is zero. For example at 
the intersection of the subspaces with GI and FD the exponents zu = z± = 3/2 
are obtained from the exponent identities. Dynamic RG recursion relations can 
be computed to one-loop order 45 ' 48 , by standard methods of momentum-shell 
dynamic RG 38 ' 39 . 

The renormalization of the seven parameters in Eqs.(52), generalized to 
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n transverse directions, give the recursion relations 



diq 

dt 

dK ± 
~df 



dt 
dT ± 
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dl 
dX ± 
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(60) 



A|| [C|| + z 
X± 



2], 



1 A,, - A 



= A> 



2C-L — C|| + Z-2 + -- 



{(X X T ± /K ± ) - (X±T ll /K ll )) 
X«K ± -X ± K\ 



■((X X T ± /K ± ) 



ir K ± K\\(K ± + K\\) 
-(A±T||/^||))] . 

The projections of the RG flows on the two parameter subspace shown in 
Fig. 8 are indicated by trajectories. They naturally satisfy the constraints im- 
posed by the non-perturbative results: the subspace of GI is closed under RG, 
while the FD condition appears as a fixed line. The RG flows, and the corre- 
sponding exponents, are different in each quadrant of Fig. 8, which implies that 
the scaling behavior is determined by the relative signs of the three nonlineari- 
ties. This was confirmed by numerical integrations 45,48 of Eqs.(52), performed 
for different sets of parameters. A summary of the computed exponents are 
given in Table I. 

The analysis of analytical and numerical results can be summarized as 
follows: 

Aj^A x > : In this region, the scaling behavior is understood best. The 
RG flows terminate on the fixed line where FD conditions apply, hence £|| = 
= 1/2. All along this line, the one loop RG exponent is z = 3/2. These 
results are consistent with the numerical simulations. The measured exponents 
rapidly converge to these values, except when A^ or A x are small. 

A x = 0: In this case the equation for ry is the KPZ equation (40), thus 
C|| = 1/2 and z\\ = 3/2. The fluctuations in ry act as a strong (multiplica- 
tive and correlated) noise on r±. The one-loop RG yields the exponents 
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z±_ = 3/2, = 0.75 for Ax > 0, while a negative Aj_ scales to suggest- 
ing z± > ay. Simulations are consistent with the RG calculations for Aj_ > 0, 
yielding (± = 0.72, surprisingly close to the one-loop RG value. For A^ < 0, 
simulations indicate zi«2 and (± w 2/3 along with the expected values for 
the longitudinal exponents. 

A^ = 0: The transverse fluctuations satisfy a simple diffusion equation 
with £j_ = 1/2 and z±_ = 2. Through the term A x (d x r±) 2 /2, these fluctuations 
act as a correlated noise 39 for the longitudinal mode. A naive application of 
the results of this reference 39 give £|| = 2/3 and z\\ = 4/3. Quite surprisingly, 
simulations indicate different behavior depending on the sign of A x . For A x < 
0, Z|| w 3/2 and Q ~ 1/2 whereas for A x > 0, longitudinal fluctuations arc 
much stronger, resulting in z\\ w 1.18 and Q| w 0.84. Actually, Q increases 
steadily with system size, suggesting a breakdown of dynamic scaling, due to 
a change of sign in A^A X . This dependence on the sign of A x may reflect the 
fundamental difference between behavior in quadrants 11 and IV of Fig. 8. 

TABLE I. Numerical estimates of the scaling exponents, for various val- 
ues of model parameters for n = 1. In all cases, = K± = 1 and 
T|| = Tx = 0.01, unless indicated otherwise. Typical error bars are ±0.05 
for C, ±0.1 for zjC^. Entries in brackets are theoretical results. Exact values 
are given in fractional form. 
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Aj^ < and A x > 0: The analysis of this region (II) is the most difficult in 
that the RG flows do not converge upon a finite fixed point and \±_ — > 0, which 
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may signal the breakdown of dynamic scaling. Simulations indicate strong 
longitudinal fluctuations that lead to instabilities in the discrete integration 
scheme, excluding the possibility of measuring the exponents reliably. 

Aj^ > and A x < 0: The projected RG flows in this quadrant (IV) 
converge to the point Aj_/A|| = 1 and A x Tj_i^|| / A|| T"|| = —I. This is actually 
not a fixed point, as K\\ and K± scale to infinity. The applicability of the CH 
transformation to this point implies zu — 3/2 and 0i = 1/2. Since A^ is finite, 
z± = Z|| = 3/2 is expected, but this docs not give any information on (±. 
Simulations indicate strong transverse fluctuations and suffer from difficulties 
similar to those in region II. 

Eqs.(52) are the simplest nonlinear, local, and dissipative equations that 
govern the fluctuations of a moving line in a random medium. They can be 
easily generalized to describe the time evolution of a manifold with arbitrary 
internal (x e R d ) and external (r e R n+1 ) dimensions, and to the motion 
of curves that are not necessarily stretched in a particular direction. Since 
the derivation only involves general symmetry arguments, the given results 
are widely applicable to a number of seemingly unrelated systems. We will 
discuss one application to drifting polymers in more detail in the next lecture, 
explicitly demonstrating the origin of the nonlinear terms starting from more 
fundamental hydrodynamic equations. A simple model of crack front propa- 
gation in three dimensions 49 also arrives at Eqs.(52), implying the self-affinc 
structure of the crack surface after the front has passed. 

2.3 Drifting Polymers 

The dynamics of polymers in fluids is of much theoretical interest and has 
been extensively studied 50 ' 51 . The combination of polymer flexibility, inter- 
actions, and hydrodynamics make a first principles approach to the problem 
quite difficult. There are, however, a number of phenomenological studies that 
describe various aspects of this problem 52 . 

One of the simplest is the Rouse model 53 : The configuration of the poly- 
mer at time t is described by a vector R(x, t), where x € [0, N] is a continuous 
variable replacing the discrete monomer index (see Fig. 9). 

Ignoring inertial effects, the relaxation of the polymer in a viscous medium 
is approximated by 

dtR(x, t) = nF(R(x, t)) = Kd*R(x, t) + n(x, t), (61) 

where /x is the mobility. The force F has a contribution from interactions 
with near neighbors that are treated as springs. Steric and other interactions 
are ignored. The effect of the medium is represented by the random forces rj 
with zero mean. The Rouse model is a linear Langevin equation that is easily 
solved. It predicts that the mean square radius of gyration, i? 2 = (|R — (R)| 2 ), 
is proportional to the polymer size N, and the largest relaxation times scale 
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Figure 9. The configuration of a polymer. 



as the fourth power of the wave number, (i.e., in dynamic light scattering 
experiments, the half width at half maximum of the scattering amplitude scales 
as the fourth power of the scattering wave vector q). These results can be 
summarized as R g <~ N" and T(q) <~ q z , where v and z are called the swelling 
and dynamic exponents, respectively 54 . Thus, for the Rouse Model, v = 1/2 
and z = 4. 

The Rouse model ignores hydrodynamic interactions mediated by the 
fluid. These effects were originally considered by Kirkwood and Riscmann 55 
and later on by Zimm 56 . The basic idea is that the motion of each monomer 
modifies the flow field at large distances. Consequently, each monomer expe- 
riences an additional velocity 

(62) 

where r xx > = R(.t) — R(x') and the final approximation is obtained by replac- 
ing the actual distance between two monomers by their average value. The 
modified equation is still linear in R and easily solved. The main result is the 
speeding up of the relaxation dynamics as the exponent z changes from 4 to 
3. Most experiments on polymer dynamics 57 indeed measure exponents close 
to 3. Rouse dynamics is still important in other circumstances, such as diffu- 
sion of a polymer in a solid matrix, stress and viscoelasticity in concentrated 



30 



polymer solutions, and is also applicable to relaxation times in Monte Carlo 
simulations. 

Since both of these models are linear, the dynamics remains invariant in 
the center of mass coordinates upon the application of a uniform external force. 
Hence the results for a drifting polymer are identical to a stationary one. This 
conclusion is in fact not correct due to the hydrodynamic interactions. For 
example, consider a rodlike conformation of the polymer with monomer length 
b where d x R a = b t a everywhere on the polymer, so that the elastic (Rouse) 
force vanishes. If a uniform force E per monomer acts on this rod, the velocity 
of the rod can be solved using Kirkwood Theory, and the result is 50 

v = t^E-p + tt]. (63) 

In the above equation, rj s is the solvent viscosity, t is the unit tangent vector, 
k = 2b/boN is the ratio of the width b to the half length boN/2 of the polymer. 
A more detailed calculation of the velocity in the more general case of an 
arbitrarily shaped slender body by Khayat and Cox 58 shows that nonlocal 
contributions to the hydrodynamic force, which depend on the whole shape 
of the polymer rather than the local orientation, are C(l/(ln k) 2 ). Therefore, 
corrections to Eq.(63) are small when N » b/b . 

Incorporating this tilt dependence of polymer mobility requires adding 
terms nonlinear in the tilt, d x r, to a local equation of motion. Since the overall 
force (or velocity) is the only vector breaking the isotropy of the fluid, the 
structure of these nonlinear terms must be identical to eq.(50). Thus in terms 
of the fluctuations parallel and perpendicular to the average drift, we again 
recover the equations, 

An A 2 

dtR\\ = U\\+K\\d 2 x R\\ + -j-(d x R\\) 2 + -^J2( d - R ^) 2 +V\\(x,t), 

i—l v ' 

d t R±i = K±dlR ±i + \±d x R\\d x R ±i + rj±i(x, t), 

where {_L i} refers to the 2 transverse coordinates of the monomer positions. 
The noise is assumed to be white and gaussian but need not be isotropic, i.e. 

f ( m (x,t) m (x',t , )) = 2T ]l S(x-x')S(t-t'), 
\ (rj± i (x,t)i] ±j (x / ,t')) =2T±SijS(x-x')S(t-t'). 

At zero average velocity, the system becomes isotropic and the equations of 
motion must coincide with the Rouse model. Therefore, {A|| , A x , Aj_, U, — 
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K±, Tii — Tj_} arc all proportional to E for small forces. The relevance of these 
nonlinear terms are determined by the dimensionless scaling variable 




where U* is a characteristic microscopic velocity associated with monomer 
motion and is roughly 10-20 m/s for polystyrene in benzene. The variable 
y is proportional to another dimensionless parameter, the Reynolds number 
i?e, which determines the breakdown of hydrodynamic equations and onset 
of turbulence. However, typically Re -C y, and the hydrodynamic equations 
are valid for moderately large y. Eqs. (64) describe the static and dynamical 
scaling properties of the nonlinear and anisotropic regime when U > U* N~ x / 2 . 

Eq.(64) is just a slight variation from (52), with two transverse components 
instead of one. Thus, the results discussed in the previous lecture apply. A 
more detailed calculation of the nonlinear terms from hydrodynamics 59 shows 
that all three nonlinearities are positive for small driving forces. In this case, 
the asymptotic scaling exponents arc isotropic, with v = 1/2 and z = 3. 
However, the fixed points of the RG transformation are in general anisotropic, 
which implies a kinetically induced form birefringence in the absence of external 
velocity gradients. This is in marked contrast with standard theories of polymer 
dynamics where a uniform driving force has essentially no effect on the internal 
modes of the polymer. 

When one of the nonlinearities approaches to zero, the swelling exponents 
may become anisotropic and the polymer elongates or compresses along the 
longitudinal direction. However, the experimental path in the parameter space 
as a function of E is not known and not all of the different scaling regimes 
correspond to actual physical situations. The scaling results found by the RG 
analysis are verified by direct integration of equations, as mentioned in the 
earlier lectures. A more detailed discussion of the analysis and results can be 
found in our earlier work . 

In constructing equations (64), we only allowed for local effects, and ig- 
nored the nonlocalities that are the hallmark of hydrodynamics. One conse- 
quence of hydrodynamic interactions is the back-flow velocity in Eq.(62) that 
can be added to the evolution equations (64). Dimensional analysis gives the 
recursion relation 

ftL =y [ vz - 1 -( d -2)v] + 0(>y 2 ), (66) 

which implies that, at the nonlinear fixed point, this additional term is sur- 
prisingly irrelevant for d > 3, and z — 3 due to the nonlinearities. For d < 3, 
z = d due to hydrodynamics, and the nonlinear terms are irrelevant. The situ- 
ation in three dimensions is unclear, but a change in the exponents is unlikely. 
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Similarly, one could consider the effect of self-avoidance by including the force 
generated by a softly repulsive contact potential 



The relevance of this term is also controlled by the scaling dimension yt — 
vz — 1 — (d — 2)f, and therefore this effect is marginal in three dimensions at 
the nonlinear fixed point, in contrast with both Rouse and Zimm models where 
self-avoidance becomes relevant below four dimensions. Unfortunately, one is 
ultimately forced to consider non-local and nonlinear terms based on similar 
grounds, and such terms are indeed relevant below four dimensions. In some 
cases, local or global arclength conservation may be an important considera- 
tion in writing down a dynamics for the system. However, a local description 
is likely to be more correct in a more complicated system with screening effects 
(motion in a gel that screens hydrodynamic interactions) where a first prin- 
ciples approach becomes even more intractable. Therefore, this model is an 
important starting point towards understanding the scaling behavior of poly- 
mers under a uniform drift, a problem with great technological importance. 
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